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INTRODUCTION 


A  multi  wavelength  laser  can  be  defined  as  a  laser  system  emitting  more  than  one 
discrete  wavelength  simultaneously.  If  each  wavelength  is  considered  a  separate  “mode” 
of  operation,  then  each  of  these  modes  may  have  differing  beam  diameters,  divergence 
values  or  even  time-intensity  variations.  In  some  instances  these  multi  wavelength 
systems  are  the  result  of  gas  lasers  operating  in  “all-line”  configurations,  or  pulsed  lasers 
equipped  w  ith  non-lmear  optical  devices  which  may  produce  fundamental,  frequency- 
doubled  and  frequency-tripled  output. 

Safety  parameters  for  multi  wavelength  lasers  cannot  be  computed  by  using  only  a  single 
analytical  mathematical  formula.  The  overall  exposure  is  increased  compared  to  that 
when  only  a  single  wavelength  is  emitted  at  a  time,  and  relative  exposure  magnitudes 
may  vary  as  a  function  of  exposure  distance  or  time.  This  prevents  applying  the  common 
equations1  for  single- wavelength  exposures  to  a  multi  wavelength  laser  problem.  While 
previous  work  by  Lyon‘  has  described  the  basic  methodology  required  to  provide  an 
approach  to  safety,  it  stops  short  of  addressing  all  considerations.  This  paper  extends  the 
methods  of  Lyon2  to  determine  optimal  laser  eye  protection  (LEP)  design,  which  both 
minimizes  overall  attenuation  and  provides  for  optimization  of  user  performance  through 
the  maximization  of  other  metrics  such  as  visible  light  transmittance. 

Single-wavelength  hazard  analysis  requires  finding  maximum  permissible  exposure 
(MPE)  limits,  accessible  emission  limits  (AEL)  and  the  effective  amount  of  energy 
passing  through  the  limiting  aperture  for  each  mode  (Qf).  This  paper  assumes  that  the 
reader  is  familiar  with  the  techniques1  for  determining  these  parameters  according  to  the 
ANSI  Z136.1-20001  Safety  Standard  and  has  a  basic  understanding  of  laser  safety 
analysis  methods  and  terminology.  Similar  constructs  and  methods  are  used  as  part  of 
the  process  outlined  here  for  multi  wavelength  analysis. 

Here,  we  examine  two  approaches  to  the  estimation  of  optical  density  (OD)  requirements 
in  a  multi  wavelength  laser  exposure.  The  first  approach  formulates  the  problem  at  a 
relatively  high  level  using  mathematical  programming,  formulated  as  a  nonlinear 
optimization  algorithm,  and  solved  using  GAMS4,  a  modeling  language  for  representing 
the  mathematical  model.  The  second  method  accomplishes  these  goals  through  the 
development  of  a  simplified  algorithm,  which  we  have  implemented  in  the  C++ 
programming  language*  for  the  purposes  of  presenting  comparative  results. 

PROBLEM  DESCRIPTION 

As  described  in  Lyon  ,  when  one  is  exposed  to  a  laser  emitting  multiple  wavelengths 
simultaneously,  the  combined  effect  should  be  considered.  In  order  to  denote  the  various 
parameters  required  for  the  analyses,  we  will  adopt  a  convention  which  is  similar  to  the 
ANSI  Z136.1  Standard.  Let  X  e  SM  (subset  of  modes)  with  the  individual  wavelengths 
emitted  by  the  laser  indexed  by  the  symbol  X.  The  symbol  Qfx  denotes  the  energy  (in 
Joules)  delivered  through  the  corresponding  limiting  aperture  at  wavelength  X  e  SM,  and 
AELx  denotes  the  accessible  emission  limit  (in  Joules)  for  wavelength  X  e  SM.  This 


1 


accessible  emission  limit  is  computed  from  the  product  of  the  MPE,  and  the  area  of  the 
appropriate  limiting  aperture.  In  absence  of  LEP,  the  safe  exposure  limit  is  considered  to 
be  exceeded  if: 


XtiSM 


(1) 


To  avoid  the  situation  described  in  Equation  (1),  LEP  is  used,  and  in  determining  optimal 
LEP  we  effectively  select  an  optical  density  Dx  at  each  wavelength  X  e  SM,  which  in 


turn  yields  a  transmittance  rx  via  tx  =  10  D*  ,X  €  SM.  To  ensure  a  safe  level  of  exposure 
wc  require: 


(2) 


Equation  (2)  will  form  one  of  the  constraints  used  in  the  mathematical  program  used  to 
optimally  design  LEP.  The  goal  is  to  provide  various  criteria  by  which  to  optimize  the 
LEP,  while  providing  safe  exposure  to  the  multi  wavelength  laser.  One  option  is  to 
maximize  visibility  through  the  eyewear,  and  for  that  we  require  a  scalar  performance 
measure  that  captures  quality  of  vision  as  a  function  of  transmittance  at  each  wavelength. 
The  spectral  luminous  efficiency  assigns  a  weight  to  each  wavelength  in  the  spectrum  that 
captures  its  relative  contribution  to  visibility  in  a  typical  environment  for  a  user*. 

There  are,  theoretically,  infinite  combinations  of  OD  for  each  multi  wavelength  laser  that 
would  provide  a  solution  to  Equation  (2).  There  may  be  some  modes  that  arc  already 
considered  safe  individually,  i.e.:  ( Qfx  / AELX )  <  1 .  For  the  other  modes,  the  OD  applied 

for  that  mode  has  a  minimum  such  that  10"Dj  ■  (< Qfx/AELx )  =  1 .  Beyond  this,  there  arc  no 
constraints  on  how  wc  apply  optical  densities. 

From  a  more  practical  standpoint,  there  may  be  other  considerations  such  as  the  cost 
involved  in  producing  OD  values  for  specific  wavelengths  or  the  desire  to  maximize 
transmittance  over  some  wavelength  band.  In  this  case,  we  may  correspondingly  assign  a 
positive  weighting  value  to  each  mode  Vx,XeA.  If  we  are  not  concerned  with 
discriminating  between  the  wavelengths  being  emitted,  a  uniform  weight  of  1.0  may  be 
applied  for  each  V x ,  which  would  in  turn  simply  minimize  the  sum  of  the  optical 
densities  for  each  mode. 
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EVAULATION  PARAMETERS 


Several  parameters  must  first  be  gathered  in  order  to  determine  the  optimal  OD 
requirement  for  each  mode  of  a  multi  wavelength  laser  exposure. 

Qfx~  The  total  amount  of  energy  expected  to  enter  the  applicable  limiting  aperture  for 

each  mode.  This  value  takes  all  modifiers  into  account  besides  OD,  including  but  not 
limited  to  atmospheric  attenuation,  transmittance  and  magnification  inherent  to  aided 
viewing  scenarios,  beam  spread  due  to  laser  geometry,  exposure  duration,  etc.  The  ANSI 
Z136.1-20001  provides  guidelines  and  calculations  for  finding  this  value  given  a 
particular  or  a  worst-case  scenario. 

AELX  -  The  accessible  emission  limit,  or  desired  energy  exposure  considered  to  make 

this  particular  mode  of  the  laser  safe.  The  ANSI  Z136.1-20001  may  be  used  to  find  the 
MPE  and  use  that  to  determine  the  worst-case  AEL  for  the  scenario.  A  less  conservative 
value  may  also  be  used  here  depending  on  the  user’s  needs. 

The  values  for  AELX  and  Qfx  may  have  been  determined  in  Watts,  Joules,  or  some  other 
related  unit  of  power  or  energy.  The  units  used  are  not  important,  so  long  as  the  units 
used  are  consistent  between  AELX  and  Of  x  such  that  the  ratio  between  them, 

( Of J  AEL ,),  is  meaningful. 

Vx  -  A  relative  weight  for  each  wavelength  that  will  help  to  discriminate  the  assignment 

of  OD  values.  This  value  must  be  greater  than  zero  at  each  wavelength,  and  must  be 
constant  with  respect  to  the  OD  assigned.  Modes  with  a  relatively  higher  V x  will  be 

assigned  a  lower  OD  such  as  to  maintain  the  objective  function.  Equation  (3).  Regardless 
of  this  weight,  the  final  solution  will  always  provide  a  safe  exposure  limit  if  the 
wavelength  were  considered  independently,  i.e.  10  Dx  ■  ( Qfx/  AELX)<  1 .  If  we  are  not 
concerned  with  discriminating  between  the  modes,  a  weight  of  1 .0  may  be  applied  for 
each  V x  ,  which  would  simply  provide  a  solution  that  minimizes  the  sum  of  the  optical 
densities  for  each  wavelength. 

An  example  weighting  value  V x ,  which  represents  the  spectral  luminous  efficiency  curve 

shown  below  in  Figure  1  as  a  function  of  wavelength,  ’  can  be  used  to  maximize 
visibility.  Because  the  human  eye  senses  different  colors  (wavelengths)  with  varying 
acuity,  when  the  wavelengths  are  outside  the  visible  spectrum,  visual  acuity  is  effectively 
zero,  but  peaks  around  550  nm.  In  our  research,  we  have  the  option  of  weighting  the 
transmittance,  t x  ,  at  each  wavelength  by  V x  and  therefore  maximize  visibility.  We  note 

that  if  these  values  are  used,  it  is  important  to  provide  some  minimal  value  for 
wavelengths  on  the  “wings”  of  the  curve  that  is  greater  than  zero.  The  examples  used  in 
the  paper  will  use  values  from  V x  shown  in  Figure  1. 
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Figure  1.  Visible  Luminous  Efficiency  for  the  Human  Eye 


PROCEDURE  FOR  OPTICAL  DENSITY  DETERMINATION 

Optimal  design  of  LEP  for  multi  wavelength  lasers  begins  with  the  same  calculations  and 
requirements  as  outlined  in  the  ANSI  Z136.1-20001  for  single  wavelength  exposures. 
However,  as  previously  indicated,  when  laser  exposures  involve  more  than  one 
wavelength,  the  calculations  are  no  longer  straightforward,  and  so.  in  this  paper 
mathematical  programming  and  algorithmic  approaches  are  used  to  determine  optimal 
LEP  design. 

Mathematical  Programming  Approach 

Mathematical  programming  is  a  set  of  models,  methods  and  theory  that  deals  with 
optimization  of  typically  complex  systems.  The  field  includes  linear  programming, 
nonlinear  programming  and  integer  programming.  It  primarily  concerns  minimizing  or 
maximizing  an  objective  function  subject  to  one  or  more  constraints.  Generally, 
mathematical  programs  seek  to  optimize  some  objective  through  the  selection  of  a  set  of 
values  for  the  decision  variables8.  These  variables  are  constrained  by  various  conditions 
and  restrictions.  Mathematical  programming  chooses  levels  of  the  decision  variables  so 
that  the  objective  is  optimized  while  maintaining  the  integrity  of  the  constraints. 

The  problem  of  optimal  LEP  design  is  first  formulated  as  a  nonlinear  optimization 
algorithm  and  solved  using  GAMS4.  An  objective  function  of  this  specific  problem 
involves  finding  a  solution  which  minimizes  OD  for  each  wavelength  and  maximizes 
visibility,  while  maintaining  the  constraint  ensuring  a  safe  exposure  scenario.  Although 
the  logarithmic  (not  linear)  relationship  between  OD  and  transmittance  appears  to  cause 
significant  difficulties,  the  specifics  of  this  problem  alleviate  this  complexity. 

The  problem  was  first  formulated  as  a  mathematical  programming  model.  Because  the 
constraint  that  ensures  a  safe  viewing  condition,  Inequality  (4),  is  nonlinear  in  form,  the 
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formulation  is  considered  a  nonlinear  program.  The  decision  variable  OD\  is  optimized 
according  to  the  following  model  formulation: 


minimize 


subject  to 


Dx  >  0,T  e  A. 


and 


(5) 


The  objective  function.  Equation  (3),  seeks  to  optimize  the  OD  according  the  user- 
specified  weighting  values  for  V When  Vx  is  the  spectral  luminous  efficiency,  the 
visibility  of  the  LEP  is  maximized  because  the  spectral  luminous  efficiency  data,  Vx, 
positively  weights  transmittance  tx  =  10'Dx  over  the  visible  wavelengths.  In  this 

situation,  the  objective  function  prefers  lower  optical  densities  in  the  visible  wavelengths 
and  therefore  more  visible  light  will  be  transmitted  through  the  LEP  to  the  user’s  eye. 

Inequality  (4)  is  the  constraint  that  ensures  the  overall  exposure  is  safe.  As  described 
previously,  a  multi  wavelength  laser  exposure  is  safe  when  Inequality  (4)  holds  for  the 
wavelengths  being  emitted  ( SM  c  A).  When  multiplied  by  the  transmittance  tx  =  1 0~Di , 

the  ratio  of  total  energy  delivered  to  the  accessible  emission  limit  must  be  less  than  one. 
It  is  through  this  constraint  that  we  have  flexibility  in  setting  the  optical  densities.  For 
example,  when  maximizing  visibility'  we  can  compensate  for  lower  optical  densities  in 
the  visible  spectrum  with  higher  optical  densities  in  the  non-visible,  or  out-of-band 
wavelengths. 

GAMS4  was  configured  using  the  formulation  above  and  run  against  a  number  of 
representative  laser  exposures.  The  results  of  the  analysis  and  their  validation  are 
described  in  Section  5. 

Algorithmic  Approach 

In  an  effort  to  utilize  the  mathematical  analysis  formulated  and  run  in  GAMS4,  an 
algorithm  was  independently  formulated  for  the  specific  problem  of  hazard  analysis. 

This  algorithm  is  sufficiently  simple  such  that  an  analysis  can  be  conducted  with  a 
handheld  calculator,  or  w'ith  basic  programming  skills  can  be  implemented  for  use  in 
software  applications.  The  algorithm  is  limited  in  scope  when  compared  to  the 
capabilities  of  GAMS,  but  should  suffice  for  most  laser  safety  analysis  requirements.  It 
is  instructive  to  examine  the  development  of  our  algorithm  through  the  use  of  an 
illustrative  example.  Working  through  the  exercise  provides  insight  to  the  logic  applied 
to  derive  the  final  equations  of  the  method. 
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After  assembling  all  of  the  output  parameters  of  the  laser  in  question,  the  initial  step  in 
the  algorithm  is  to  determine  if  the  laser  is  already  safe  without  any  additional  OD  from 
LEP.  This  involves  simply  summing  up  the  ratios  of  each  ( Qfx/  AELx)and  checking  to 

see  if  the  sum  is  less  than  or  equal  to  one,  i.e.  ^  \Qfx  / AELX )  <  1 .  The  remainder  of  the 

x 

procedure  assumes  that  this  is  not  the  case.  The  goal  of  this  algorithm  then  is  to  fulfill  (3) 
and  (4),  restated  as: 

A  =  1  0  A  ^ —  =  ]  whjie  minimizing  Z  =  ^  '  D x  •  Vx . 

ZELX  x 

If  we  consider  these  equations  for  each  mode  (wavelength)  independently  and  take  the 
first  derivative  of  each  with  respect  to  Dx,  we  can  easily  see  the  effect  of  increasing  the 
OD  for  each  mode: 


^=10 


-D»  Qfx 


AEL- 


A'  =1(Ta  - ln(10) 

AELX 


(6) 


=  Dx  ■  Vx 

Z'x=Vx.  (7) 


Therefore,  our  approach  is  to  set  the  goal  of  the  algorithm  to  increase  OD  on  the  mode 
that  has  the  greatest  impact  in  decreasing  A  to  1  while  having  the  least  impact  in  raising 
Z.  To  help  us  do  that,  we  will  introduce  a  weighting  function  that  combines  (6)  and  (7): 


Wx(Dx)  =  -^  =  10-* 


Qfx 

AELX  •  Vx 


-ln(10). 


(8) 


Since  we  are  only  using  the  function  for  comparison  between  modes,  wc  can  factor  out 
the  common  constant  elements  to  get: 


WX(DX)  =  10  A- 


Qfx 

AELx  •  Vx  ‘ 


(9) 


The  most  beneficial  effect  of  adding  small  increments  of  OD  to  individual  modes  of  the 
laser  is  adding  OD  to  the  mode  with  the  greatest  IT  value  for  its  current  Dx.  To  help 
illustrate  this,  we  will  use  the  following  example  of  a  multi-mode  laser  (this  is  the  laser 
named  “C”  in  the  results  section).  We  have  also  used  the  spectral  luminous  efficiency 
curve  described  above  to  obtain  illustrative  values  for  V.  Table  1  illustrates  the  initial 
values  for  the  algorithm  using  laser  “C”. 
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Table  1.  Initial 

Algorithm  Values  for  Example  Laser  “C” 

Of 

AEL 

V 

Wavelength 

(W) 

(W) 

(nm) 

458 

0.02 

6.50E-04 

5.56E-02 

488 

0.1 

6.50E-04 

1.94E-01 

514 

0.2 

6.50E-04 

5.86E-01 

640 

0.05 

6.50E-04 

1.75E-01 

Examination  of  Table  1  quickly  shows  that  each  of  the  modes  of  this  laser  is  hazardous 
(Qf/AEL  >  1  for  each  mode),  so  each  mode  has  a  minimal  OD  value  ( Dmin )  greater  than 
zero.  The  values  of  Dmin  for  each  mode  are  shown  in  Table  2. 


Table  2.  Individual  Mode  Hazard  Analysis  Values  of  Dmin  and  Weighting  Factors  for  Example  Laser  “C” 


Wavelength 

Qf 

AEL 

V 

T 

Dmin 

W 

(nm) 

(W) 

(W) 

458 

0.02 

6.50E-04 

5.56E-02 

3.25E-02 

1.49 

18.0 

488 

0.10 

6.50E-04 

1.94E-01 

6.50E-03 

2.19 

5.15 

514 

0.20 

6.50E-04 

5.86E-01 

3.25E-03 

2.49 

1.71 

640 

0.05 

6.50E-04 

1. 75E-01 

1.30E-02 

1.89 

5.71 

Applying  a  minimal  OD  of  Dx  from  Table  2  makes  each  mode  safe  (Ax  for  each  mode  is 

1),  but  the  laser  as  a  whole  is  still  not  safe  because  there  are  four  modes.  We  see  that 
Equation  (6)  would  yield  a  value  of  4  in  this  case.  We  still  need  to  increase  the  OD  for 
one  or  more  modes.  We  can  choose  the  mode  to  work  on  by  computing  W(D). 

If  we  solve  for  W(Dmin)  for  each  mode,  we  can  find  the  mode  with  the  greatest  IE  that  will 
allow  us  to  increase  its  OD  with  the  greatest  impact  on  A  and  the  least  impact  on  Z. 

Table  2  illustrates  the  computation  of  W  for  each  mode  based  upon  the  initial  hazard 
analysis. 

From  the  values  in  Table  2,  we  see  that  the  mode  to  address  is  the  one  with  the  highest  W 
value  based  on  Dmm.  Raising  the  OD  for  the  458-nm  mode  a  small  amount  will  have  the 
greatest  impact  on  A  and  the  least  impact  on  T.  But  how  much  can  we  raise  this  OD? 
Looking  at  Equation  (9),  we  will  note  that  W  is  monotonically  decreasing  as  D  increases 
for  all  values  that  are  greater  than  or  equal  to  0.  This  means  continually  raising  the  OD 
for  one  mode  will  have  a  decreasingly  beneficial  effect.  We  can  increase  the  OD  of  the 
458-nm  mode  until  W(D)  is  equal  to  the  next  lowest  IE  value;  in  this  case,  the  640-nm 
mode  has  a  IE  value  of  5.71. 
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Using  Equation  (9)  and  solving  for  Dx  we  obtain  Equation  (10). 


( 

Da=~  log 

V 


K  •  ael,  •  vK 
Qf. \ 


/ 


(10) 


In  this  case,  we  get  a  new  value  of  for  the  458-nm  mode  of  1 .99.  We  now  compute 
new  values  of  A  and  Z  to  determine  if  the  new  combination  of  OD  values  will  make  the 
entire  laser  safe.  Table  3  represents  an  updated  configuration  of  parameters,  along  with 
the  values  of  A  and  Z. 


Table  3.  Individual  Mode  Hazard  Analysis  Values  of  Dmin  and  Weighting  Factors  for 
Example  Laser  “C”  after  One  Algorithm  Iteration 


Wavelength 

(nm) 

Qf 

(W) 

AEL 

(W) 

V 

Dx 

W 

458 

0.02 

6.50E-04 

5.56E-02 

1.99 

5. 71 

488 

0.1 

6.50E-04 

1.94E-01 

2.19 

5.15 

514 

0.2 

6.50E-04 

5.86E-01 

2.49 

1.71 

640 

A  =  3.3177,  Z  = 

0.05 

2.3288 

6.50E-04 

1. 75E-01 

1.89 

5.71 

Unfortunately,  this  still  does  not  yet  provide  for  a  safe  exposure  (A  >  1),  so  we  need  to 
continue.  We  can  now  look  to  the  next  highest  If 'value,  and  increase  the  OD  at  the  458- 
nm  and  the  640-nm  modes  to  the  point  where  their  IE  values  are  the  same.  In  this  ease, 
the  488-nm  mode  has  a  ff'value  of  5.15.  We  ean  apply  (10)  to  the  458-nm  and  540-nm 
modes  to  find  new  Dx  values  for  both  of  them.  The  results  are  shown  below'  in  Table  4. 
We  see  that  the  value  of  A  is  indeed  reduced,  indicating  that  the  exposure  is  less 
hazardous.  However,  A  remains  greater  than  1 ,  indicating  that  a  hazard  does  indeed 
remain. 


Table  4.  Individual  Mode  Hazard  Analysis  Values  of  Dmin  and  Weighting  Factors  for 
Example  Laser  “C”  after  Two  Algorithm  Iterations 


Wavelength 

(nm) 

Qf 

(W) 

AEL 

(W) 

V 

Dx 

W 

458 

0.02 

6.50E-04 

5.56E-02 

2.03 

5.15 

488 

0.1 

6.50E-04 

1.94E-01 

2.19 

5.15 

514 

0.2 

6.50E-04 

5.86E-01 

2.49 

1.71 

640 

A  =  3. 1874, Z 

0.05 

=  2.3332 

6.50E-04 

1.75E-01 

1.93 

5.15 
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At  every  stage  of  iteration,  we  are  using  W  to  find  the  optimal  mode  and  only  increasing 
the  OD  for  that  mode  which  has  the  greatest  impact  in  decreasing  A  with  the  least  impact 
on  increasing  Z.  We  could  continue  with  this  method  until  we  reach  A  =  1. 


We  see  that  if  we  had  to  increase  the  OD  for  every  mode,  then  every  mode  would 
eventually  share  the  same  W value.  This  indicates  that  we  can  solve  for  this  W value  (W) 
from  the  very  beginning  and  use  Equation  (10)  to  find  the  OD  that  would  provide  that  W’ 

Of 

value.  Using  Z  =  Y'lO"£>'1 - —  =  1  and  replacing  OD  using  (10)  we  get: 

J  4  £  T 


5> 


A 


Qfx  _x 

AELX 


(ii) 


This  equation  simplifies  to: 

ay 

X  X 

for  cases  of  constant  values  of  W\  This  indicates  that  we  can  again  simplify  to: 


2A 

A 

(13) 


So  for  our  example  laser  “C”,  we  can  find  W 9  =  0.99  from  our  values  of  Vx.  This  gives 
the  resulting  optical  densities  shown  below  in  Table  5,  and  a  safe  multi  wavelength  laser 
exposure. 


Table  5.  Individual  Mode  Hazard  Analysis  Values  and  Weighting  Factors  for 
Example  Laser  “C”  Applying  Equation  (13)  in  the  Development  of  our  Algorithm 


Wavelength 

Qf 

AEL 

V 

(nm) 

(W) 

(W) 

458 

0.02 

6.50E-04 

5.56E-02 

488 

0.1 

6.50E-04 

1.94E-01 

514 

0.2 

6.50E-04 

5.86E-01 

640 

0.05 

6.50E-04 

1.75E-01 

A  =  1.0,  Z  = 

2.7762 

Dx 

W 

2.75 

0.99 

2.90 

0.99 

2.72 

0.99 

2.65 

0.99 
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For  a  multi-mode  laser  where  every  mode  is  hazardous,  the  solution  then  is  to  find  W’ 
using  Equation  (13),  then  find  the  value  oiDx  for  each  mode  using  Equation  (10).  If 
there  are  one  or  more  modes  that  are  not  hazardous,  we  may  not  know  if  those  modes 
require  any  OD  as  part  of  the  safety  analysis.  As  an  example  of  this  condition,  we  take 
the  example  laser  “C”  and  change  the  value  of  Qf  for  the  640-nm  mode  to  be  non- 
hazardous  and  evaluate  the  optimal  OD  requirements.  We  begin  with  the  parameters  in 
Table  6.  We  can  assume  an  initial  OD  of  zero  for  each  mode  and  use  Equation  (9)  to 
solve  for  the  W  values  for  each  mode.  This  analysis  is  also  shown  in  Table  6. 


Table  6.  Initial  Parameters  for  a  Modified  Laser  “C”  Assessment  of  Optical  Density  Requirements 


Wavelength 

Qf 

AEL 

(nm) 

(W) 

(W) 

458 

0.02 

6.50E-04 

488 

0.1 

6.50E-04 

514 

0.2 

6.50E-04 

640 

6.0E-04 

6.50E-04 

V 

Dmin 

w 

5.56E-02 

0 

553.4 

1.94E-01 

0 

793.0 

5.86E-01 

0 

525.1 

1. 75E-01 

0 

5.3 

In  this  situation,  we  will  note  that  the  W(0)  values  for  the  mode  are  still  greater  than  the 
W’  value  calculated  using  Equation  (13),  which  was  0.99.  This  means  that  even  though 
the  640-nm  mode  is  not  hazardous,  the  optimal  solution  would  still  involve  raising  the 
OD  of  this  mode  some  amount.  Once  again,  we  apply  Equation  (10)  to  calculate  the  OD 
for  each  mode  in  order  to  provide  for  a  safe  exposure.  We  obtained  the  results 
summarized  by  Table  7. 


Table  7:  Algorithm  Results  for  Assessment  of  a  Modified  Laser  “C” 


Wavelength 

Qf 

AEL 

V 

Dx 

W 

(nm) 

(W) 

(W) 

458 

0.02 

6.50E-04 

5.56E-02 

2.75 

0.99 

488 

0.1 

6.50E-04 

1.94E-01 

2.90 

0.99 

514 

0.2 

6.50E-04 

5.86E-01 

2.72 

0.99 

640 

6.00E-04 

6.50E-04 

1. 75E-01 

0.73 

0.99 

A  =  1.0,  Z=  2.4397 


If  one  or  more  of  the  W(0)  values  for  the  mode  are  less  than  or  equal  to  the  W’  value,  then 
these  modes  do  not  need  any  OD  to  reach  an  optimal  solution.  When  this  occurs,  it  is 
appropriate  to  remove  those  modes  from  the  equation  and  solve  for  a  new  W’  value.  This 
can  be  done  until  all  of  the  remaining  modes  have  a  JV(O)  value  of  greater  than  W’.  As  a 
final  case  study,  we  modify  our  example  laser  “C”  to  present  the  situation  when  no  OD  is 
required  at  one  wavelength.  Our  initial  parameters  are  presented  in  Table  8. 
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Table  8.  Initial  Parameters  and  for  a  Second  Modified  Laser  “C”  Assessment 
of  Optical  Density  Requirements 


Wavelength 

(nm) 

Qf 

(W) 

AEL 

(W) 

V 

Dmin 

W 

458 

0.02 

6.50E-04 

5.56E-02 

2.71 

1.08627 

488 

0.1 

6.50E-04 

1.94E-01 

2.86 

1.08627 

514 

0.2 

6.50E-04 

5.86E-01 

2.68 

1.08627 

640  6.00E-05 

A  =  1.0,  Z=  2.2760 

6.50E-04 

1.75E-01 

0 

N/A 

In  this  example,  the  Qf  value  for  the  640-nm  mode  is  lowered  such  that  its  W  value  is  less 
than  the  0.99  that  was  calculated  for  W’.  This  indicates  that  the  optimal  solution  would 
not  include  any  OD  for  this  mode.  We  also  cannot  have  a  negative  OD  value,  so  we 
should  effectively  remove  this  mode  from  the  computation. 


First,  we  calculate  the  effect  that  this  mode  would  have  on  A.  In  this  case, 

A ,  =  1 0~D 1  ■  Qfx  / AEL ,  for  the  640-nm  mode,  such  that  Ax  -  — — — — — —  =  0.0923 1 .  We 

A  x  6.50E-4 

effectively  want  the  sum  of  Ax  of  the  other  three  modes  to  make  the  difference  between 

1.0  and  0.0923,  or  0.9077.  Using  this  new'  value  for  our  desired  A  ’,  w'e  will  find  the  new 

optimal  W'  using  a  derivation  of  Equation  (13):  W  =  0.90769/  Y/x  for  all  of  the 

/  /UA 

remaining  modes.  So  in  this  case  we  have:  W'  =  1 .08627 .  We  then  use  Equation  (10)  to 
find  the  OD  for  each  of  the  remaining  modes,  giving  us  the  optimal  solution  shown  in 
Table  9. 


To  re-cap,  the  procedure  can  be  summarized  as  follows: 


(1)  Find  Wx  for  each  mode  of  the  laser,  starting  with  an  OD  of  0  for  each  mode.  Using 
the  simplified  Equation  (9). 


Wx(  0)  = 


QL 

AELX-VX 


(14) 


(2)  Find  the  optimal  W’  for  the  laser  using  (13): 


ZX 


(15) 
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Compare  W  to  the  Wx  for  each  mode.  If  there  are  any  W x  <  W  ' ,  these  modes  will  not 
have  any  increase  in  Dx-  Call  this  subset  of  A,  J  and  remove  them  from  the  calculation 
of  W  using: 


V . 

jLu 


w  = 


r  aelj 


(16) 


If  this  new  W’  is  still  greater  than  Wx  of  one  of  the  remaining  modes,  remove  that  mode 
in  the  same  way  as  the  step  above  until  there  is  some  sub-set  of  the  original  modes  where 
Wx  >  W  for  each  mode. 


(3)  For  every  mode  remaining,  find  the  required  optical  density  using: 


( 


Dx  =-log 


W  ■  AELX 

~qK 


•  V, 


\ 

J 


(17) 


We  note  that  in  case  of  equally-weighted  10. =1.0,  then  the  trivial  solution  of  increasing 
each  optical  density  by  logio(AO,  where  A  is  the  number  of  individual  laser  modes. 
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DISCUSSION 


The  two  methods  were  checked  for  agreement  and  accuracy  using  data  from  five  example 
lasers;  the  results  are  shown  in  Table  9.  The  first  column  of  the  table  below  contains  the 
wavelengths  that  each  example  laser  emits.  The  second  and  third  columns  are  the  Qfx 

and  AELX  as  calculated  according  to  the  ANSI  Z136.1-20001.  The  resulting  optimal 

optical  densities  at  each  wavelength  are  shown  for  the  two  methods:  nonlinear 
programming  and  the  algorithmic  approach  in  the  fourth  and  fifth  columns.  The  safety 
check  is  a  manual  calculation  of  the  value  of  the  constraint  required  for  a  safe  exposure 
from  Inequality  (2).  Recall  that  the  value  of  the  safety  constraint,  Inequality  (2),  must  be 
less  than  or  equal  to  1  in  order  to  ensure  a  safe  exposure. 


Table  9.  Summary  of  Data  Comparing  Mathematical  Programming  Results 
and  the  Algorithm  for  Five  Sample  Lasers 


Laser  Input  Paramters 


Optical  Density  Results 


“A” 

Qf(J) 

532  nm 

1.00E-03 

1 064  nm 

5.00E-04 

“B” 

Qf(J) 

532  nm 

1.00E-01 

860  nm 

5.00E-01 

“C” 

Qf(J) 

458  nm 

2.00E-02 

488  nm 

1.00E-01 

514  nm 

2.00E-01 

640  nm 

5.00E-02 

“D” 

Qf(J) 

430  nm 

1.00E-02 

530  nm 

5.00E-04 

640  nm 

5.00E-02 

“E” 

Qf(J) 

800  nm 

1.00E-02 

1064  nm 

1 .00E-03 

AEL(J) 

5.00E-07 

5.00E-06 

safety  check 

AEL(J) 

6.50E-04 
1 .40E-03 

safety  check 

AEL(J) 

6.50E-04 

6.50E-04 

6.50E-04 

6.50E-04 

safety  check 

AEL(J) 

1.00E-05 

6.50E-04 

6.50E-04 

safety  check 

AEL(J) 

1.60E-02 

5.00E-06 

safety  check 


Nonlinear 

Algorithmic 

rogramming 

approach 

3.301 

3.302 

10.000 

4.945 

1.000 

0.999 

2.187 

2.188 

10.000 

5.498 

1.000 

0.999 

2.748 

2.748 

2.903 

2.903 

2.725 

2.725 

2.648 

2.648 

1.000 

1.000 

4.829 

4.829 

0.000 

0.000 

2.552 

2.552 

1.000 

1.000 

2.085 

0.097 

2.303 

2.602 

1.000 

1.000 
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The  first  results  column  presented  in  Table  9  are  based  on  solving  the  nonlinear 
program  using  the  MINOS  NLP  solver,  which  we  call  from  the  GAMS  mathematical 
programming  language9  4.  These  are  compared  to  the  results  from  the  algorithmic 
approach  implemented  in  the  C++  programming  language5,10.  We  see  by  the  overly 
safe  OD  of  1 0  that  resulted  in  a  few  of  the  examples  (with  the  nonlinear  mathematical 
program)  that  this  particular  formulation  is  not  the  best  for  that  method,  but  provides 
a  basis  and  comparison  for  the  algorithmic  approach. 

By  comparing  the  optical  densities  determined  to  be  optimal  for  each  example  laser,  it 
is  clear  that  the  wavelengths  in  the  visible  spectrum  are  being  favored  according  to 
the  weights  set  by  the  spectral  luminous  efficiency  values.  These  weights  prefer  lower 
optical  densities  for  the  visible  wavelengths  and  the  models  compensate  by  putting 
higher  OD  on  the  non-visible  wavelengths.  This  satisfies  the  goal  of  maximizing 
visibility  while  providing  a  safe  viewing  condition  for  multi  wavelength  exposures. 
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CONCLUSIONS 


This  paper  presents  a  brief  tutorial  regarding  the  computation  of  optimal  optical  density 
values  for  multi  wavelength  lasers.  It  is  our  hope  that  with  this  information,  hazard 
analyses  may  be  conducted  with  more  understanding  and  greater  confidence  in  the 
results.  The  purpose  of  this  work  was  to  transition  the  optimization  methods  used  in 
nonlinear  programming  into  a  versatile  format  that  can  be  integrated  into  existing  laser 
hazard  software,  LTMC. 

The  Laser  Threat  Modeling  Components  (LTMC)  is  an  ANSI  C++  library  developed  by 
AFRL/HEDO  to  encapsulate  the  computational  logic  used  by  LHAZ  and  LRMS  to 
conduct  basic  ANSI  Z136. 1-2000  analyses.  It  also  contains  additional  logic  for 
computing  ED50's,  hazard  distances,  recommended  optical  densities,  and 
irradiance/radiant  exposure  at  range.  It  has  been  compiled  under  Windows,  Mac  OS  X, 
and  Linux.  The  current  implementation  supports  single  wavelength  emitting  lasers.  For 
these,  as  indicated  earlier,  a  recommended  optical  density  is  computed  to  ensure  exposure 
at  a  specified  minimal  range  does  not  exceed  the  ANSI  Z136. 1-2000  maximum 
permissible  exposure  (MPE)  threshold.  Plans  for  the  next  generation  LHAZ  application 
include  support  for  multi  wavelength  laser  hazard  analyses.  In  conjunction  with  this 
plan,  an  implementation  of  the  OD  optimization  algorithms  discussed  here  will  be 
completed  as  an  extension  to  LTMC's  support.  Additionally,  a  graphically  user  interface 
within  LHAZ  will  be  supplied  to  access  this  LTMC  capability. 
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APPENDIX 


This  appendix  contains  the  C++  source  code  for  the  algorithmic  approached  detailed 
above. 


#include  "SMinOD.h" 

#include  <vector> 

#include  <math.h> 
using  namespace  std; 

III  Looks  through  a  vector  of  integers  to  see  if  index  exists 
bool  Contains (cons t  vector<int>&  UsedI,  int  index) 

{ 

for  (unsigned  i  =  0;  i  <  UsedI . si ze  () ;  ++i) 

{ 

if  (UsedI [i]  ==  index)  return  true; 

} 

return  false; 

} 

III  The  function  that  does  the  work.  A  vector  of  SMinOD  instances  is  passed 
in  . 

Ill  The  OD  member  for  each  instance  is  populated  based  on  the  algorithm 

described  in  the  documentation 

void  SMinOD :: MinimizeOD (vector<SMinOD>&  V) 

{ 

II  Watch  for  pedantic  case  of  already  being  under,  this  also  checks  for 
0  elements 

double  A1  =  1.0; 

double  safety  *  0.0; 

int  c  =  (int) V. size () ;  int  i  =  0; 

for  (i  =  0;  i  <  c;  ++i) 

{ 

V [ 0 ] . OD  =  0;  //  Populate  the  required  OD  at  0 

for  the  pedantic  case 

safety  +=  V[i].Qf  /  V[i].AEL; 

} 

if  (safety  <=  1.0)  return; 

//  Look  for  a  single  mode,  find  the  simple  OD  for  that  mode 
if  (c  ==  1) 

{ 

V [ 0 ]  . OD  =  -1.0  *  loglO  (V[0]  .AEL/V[0]  .Qf)  ; 
return; 

} 

//  Watch  for  modes  that  have  maximum  OD,  see  if  it  will  drop  the  value 
below  1 

for  (i  =  0;  i  <  c;  ++i) 

( 

if  (V [i] .VisC  >1.0) 

{ 

V [ i ] .OD  =  V [ i ] .VisC; 

safety  -=  V[i].Qf  /  V[i].AEL; 
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} 


safety  +=  (V[i].Qf  *  pow{10.0,  -V[i].OD}}  /  V[i] .AEL; 


if  {safety  >  1.0} 

{  //  Safety  is  still  high,  max  the  OD  on  what  we  have  and  use 

W  for  the  other  modes 

for  {i  =0;  i  <  c;  ++i} 

{ 

if  ( V [ i ] .OD  ==  0.0} 

V [ i ] .OD  =  -V [ i ] . Qf  /  {V[i).AEL  *  V[i].VisC}; 

else 

A1  -=  (V [ i ] . Qf  *  pow {10.0,  -V [ i ) . OD} }  / 

V [ i ] .AEL; 

} 


} 

else  if  {safety  <  1.0} 

{  //  Safety  is  TOO  low.  Make  all  of  the  other  numbers  0  and 

just  solve  for  these  maxOD's 

if  ( V [ i ] .OD  ==  0.0} 

A1  -=  V [ i ] . Qf  /  V[i] .AEL; 
else  //  Weight  based  on  the  max  OD 
{ 

V[i) .VisC  =  1.0  /  V [ i ] . OD ; 

V [ i ] . OD  =  -V [ i ] . Qf  /  (V [ i ] . AEL  *  V[i].VisC); 

} 

} 

else  return;  //  Probably  won't  happen,  but  just  in  case  max 

OD  makes  us  perfectly  safe... 


//  Loop  until  we  find  a  W  that  works  for  the  modes,  and  use  0  for  the 
OD  of  the  others 

double  W; 

double  minW  =  le6; 
do 
{ 

•  minW  =  le6; 

W  =  0.0; 

int  mini  =  -1; 

for  (i  =  0;  i  <  c;  ++i) 

{ 

if  (V [i] .OD  <  0.0) 

{ 

W  +=  V [ i ) .VisC; 
if  ( —V [ i ] .OD  <  minW} 


mini  =  i; 
minW  =  -V [ i ] .OD; 

} 

} 

} 

W  =  A1  /  W; 
if  (W  >  minW) 

{ 

V [mini ] .OD  =  0.0; 

A1  -=  V [mini ] . Qf  /  V[minl] .AEL; 

} 
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) 

while  (W  >  minW) ; 


//  Find  the  OD  for  the  other  modes 
for  (i  =  0;  i  <  c;  ++i) 


V[i] .Qf) ; 


if  ( V [ i ] .OD  <  0.0) 

V [i ]  .OD  =  -1.0  * 


) 


} 


loglO  (  (W 


V [ i ] .AEL 


V [ i ] . VisC)  / 
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